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Bimetallic rod-shaped nanomotors swim autonomously in hydrogen peroxide solutions. Here we 
present a scaling analysis, computational simulations, and experimental data that show that the 
nanomotor locomotion is driven by fluid slip around the nanomotor surface due to electrical body 
forces. The body forces are generated by a coupling of charge density and electric fields induced 
by electrochemical reactions occurring on the nanomotor surface. We describe the dependence of 
nanomotor motion on the nanomotor surface potential and reaction-driven flux. 
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Autonomously propelled synthetic nanomotors rep- 
resent a major step in the development of practical 
nanomachines because they are able to perform nanoscale 
tasks without the need for externally supplied energy. 
Synthetic nanomotors often take the form of Janus 
nanoparticles consisting of two dissimilar segments sus- 
pended in a liquid fuel with asymmetric reactions occur- 
ring on the nanomotor surface. Examples include the 
catalytic p]-[3| and electrocatalytic 0, [H| decomposition 
of hydrogen peroxide and enzymatic reactions 0,0] • The 
nature of the reaction and the locomotion mechanism at 
work depend on the fuel and particle materials. Several 
particle locomotion mechanisms have been investigated 
including autonomous variants of diffusiophoresis [6| , os- 
miophoresis_p, [1| , electrophoresis 0, 0| , surface-tension 
gradients (IP) , and bubble propulsion [2|. In addition, 
Golestanian and coworkers have provided a general math- 
ematical framework for phoretic swimmers based on a 
particle's surface activity and mobility that can be ap- 
plied to several of the aforementioned mechanisms [11] . 

Several groups have fabricated self-propelled bimetal- 
lic nanowire motors that swim due to the electrochem- 
ical decomposition of hydrogen peroxide (H2O2) fuel to 
oxygen and water 0-0 • These electrochemically grown 
nanowires (or nanorods) have been engineered to: 
swim at one hundred body lengths per second 
(ii) perform controlled motion under applied magnetic 
fields d, [l3(| , chemical [l4j and thermal [HI modula- 
tion; (iii) sense chemicals through their autonomous mo- 
tion [16]; and (iv) pick up, haul, and release micron-scale 
cargo [4| . Several of the aforementioned studies have con- 
tributed conceptual models and experimental data in an 
effort to elucidate the mechanism by which these bimetal- 
lic nanorods convert chemical energy into motion. How- 
ever, thus far there is no universally accepted theory that 
fully describes the locomotion physics. 

In this Letter, we present a set of governing equa- 
tions, a scaling analysis, numerical simulations, and ex- 
periments that describe the physics underlying the au- 
tonomous motion of electrocatalytic bimetallic nanomo- 
tors due to a mechanism we call Reaction Induced Charge 
Auto-Electrophoresis (RICA). To our knowledge, this 
is the first work which solves the Poisson, advection- 
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FIG. 1: Schematic of Pt/Au catalytic nanomotor showing 
electrochemical reactions, generated charge density (red and 
blue denote high and low charge density, respectively), and 
approximate electric field lines. Although a gold-platinum 
nanomotor is shown here, the results in this work are ap- 
plicable to other bimetallic combinations that undergo elec- 
trochemical reactions As shown here, this nanomotor's 
autonomous motion would be directed to the right. 



diffusion, and Navier-Stokes equations to provide a de- 
tailed physical description of the locomotion of bimetallic 
nanorods. The scaling analysis and simulations enable 
predictions of nanomotor velocity, direction, and total 
propulsive force as a function of two relevant parameters 
of the system. We compare these predictions to experi- 
mental measurements and observe excellent agreement. 

Figure 1 shows a schematic of a Pt/Au bimetallic 
nanomotor in an aqueous hydrogen peroxide solution. 
Peroxide oxidation at the Pt end generates protons, elec- 
trons, and oxygen molecules. The electrons conduct 
through the rod to the Au end and complete the reduc- 
tion reaction by combining with protons, peroxide and 
oxygen to generate water. This rod effectively acts like a 
short-circuited galvanic cell which drives a net migration 
of protons in the surrounding solution from the anode 
(Pt) to the cathode (Au). Wang et al. used this basic de- 
scription along with electrochemical measurements ^ to 
successfully predict the direction of motion of a nanorod 
composed of any 2 of 6 noble metals. Although there is 
a growing consensus that this mechanism is operative in 
causing the rods' motion, the underlying physical details 
are not well understood. 

Here we provide a model that shows that the asymmet- 
ric reactions result in an excess and depletion of protons 
in the surrounding electrolyte at the anode and cath- 
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ode ends, respectively. The proton imbalance results in 
asymmetric free charge density, as shown in Fig. 1, which 
generates an electric dipole and field pointing from the 
anode to the cathode. The self-generated electric field 
couples with the charge density to produce an electri- 
cal body force that drives fluid from the anode to the 
cathode. The fluid motion results in locomotion of the 
nanowire in the direction of the anode. This physics is 
similar to electrophoresis, except here the electric field 
and charge density are induced by the surface reactions. 

In our model, we consider a rod immersed in a binary 
electrolyte with equal concentrations of H + and OH - 
ions. In the dilute solution limit, the steady ion concen- 
tration distributions are given by the advection-diffusion 
equation, 
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where u is the fluid velocity, c is the molar concentration, 
D is the diffusivity, z is the valence, F is the Faraday 
constant, v is the mobility, <fi is the electrostatic potential, 
and the subscript i denotes the species. The electrostatic 
potential in turn depends on the local free charge density 
as described by the Poisson equation, 



eV 2 



F(z + c + + Z-C-) 



(2) 



where z+ = 1, Z- = —1, p e is the volumetric charge 
density, and e = e r eo is the permittivity of the liquid 
which we assume to be constant. In our model, fluid 
motion (and thus nanomotor motion) is driven by elec- 
trical body forces, which depend on the concentrations of 
charged species only. We have estimated the forces due 
to diffusiophoresis and found that they are several orders 
of magnitude smaller than those produced by the induced 
charge mechanism described here. For this reason we do 
not consider the oxygen and peroxide concentrations and 
focus on a simple binary electrolyte. To close the sys- 
tem of equations we include the steady, incompressible 
Navier-Stokes equations for a Newtonian fluid, 
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Here p is the fluid density, p is the pressure, r] is the 
dynamic viscosity, and p e V(/> is the electrical body force 
that results from the coupling of charge density and elec- 
tric field. This general framework was developed by 
Melcher and has been used extensively to describe elec- 
trohydrodynamic flows, particularly electrokinetic flows 
which are driven by a coupling of externally applied fields 
and charged objects [l7| . 

The reactions are represented by boundary conditions 
specifying the molar proton fluxes on the surface of the 
nanomotor. On the anode and cathode we prescribe 
equal and opposite fluxes j and — j normal to the wire 
surface. Since anions (hydroxide ions) do not participate 
in the reactions, the normal anion flux is set to zero ev- 
erywhere on the nanomotor surface. The values of the 



proton fluxes specified in the simulations are based on 
previously published measurements of current density at 
Pt and Au electrodes in hydrogen peroxide [l2|, [3, [HI . 
Here we do not directly model the electrochemical reac- 
tions that are described by the Butler- Volmer equation, 
because we have direct measurements of the current den- 
sity for our electrocatalysts and fuel. At the nanomotor 
surface we apply the no-slip condition for the velocity 
and specify the local surface potential (relative to the 
bulk solution) as (j) = £. Far from the rod surface, the 
electrostatic potential decays to zero and ion concentra- 
tions approach their bulk value, i.e., <fi —> and q —> 
as the radial distance r — >> oo. 

We non-dimensionalize the momentum equation using 
the following scaling quantities: \u\ oc U ev , p oc r]U ev /d, 
p e oc p eQ1 and V(/> oc Eq where p eo is a characteristic 
charge density, Eq is a characteristic electric field, d is a 
viscous length scale, and U ev is a characteristic electrovis- 
cous velocity. Applying these scalings to the momentum 
equation (4), a Reynolds number emerges based on the 
electroviscous velocity, given by Re = pU ev d/rj. Here we 
use the electroviscous velocity which arises in electroki- 
netic systems due to the balance of electrical body and 
viscous forces acting on the fluid, defined as [ItJ 
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Metallic nanorods support a native surface charge in 
aqueous solutions [19|. The charged surface attracts 
a screening cloud of counter-ions which develops a re- 
gion of net charge density in the electrical double layer 
(EDL) surrounding the rod. The characteristic length 
scale of the charge density region is the Debye thickness 
A d , which scales with the background electrolyte concen- 
tration Coo -1 / 2 . This charge density in the EDL scales 
with the potential in the EDL based on the Poisson- 
Boltzmann equation for a symmetric binary electrolyte, 
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where we have imposed the Debye-Hiickel approximation. 
In order to also include the effects of the reaction-driven 
flux, we introduce a characteristic electric field based on 
the flux and diffusivity of protons, 
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where h is the length of the nanomotor and also the char- 
acteristic length for the electric field. Combining expres- 
sions (5-7), the electroviscous velocity scales as 
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Here we scale the viscous length scale d with the EDL 
thickness since the region of significant viscous and 
electrical body forces is limited to the area with charge 
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FIG. 2: Simulation-generated plots of (left) normalized proton concentration c* (color) and electrical potential 0* (contour 
lines) and (middle) charge density p e /pe (color) and electric field (streamlines) and (right) RICA velocity magnitude (colors) 
and streamlines (black lines) for the case £ = — 10 mV and j/jd = 0.8. The reactions lead to an asymmetry in the proton 
concentration such that an excess of protons builds up at the anode and protons are depleted at the cathode. The excess 
of protons results in positive charge density at the anode and the generation of electric field pointing from the anode to the 
cathode. 



density. As will be shown below, the ^-component of 
the self-generated electric field is significant along a dis- 
tance that scales with the nanomotor length h. The scal- 
ing analysis shows that the nanomotor speed increases 
linearly with the reaction flux j because the flux gener- 
ates charge density which produces the internally gener- 
ated field. Equation (8) can be recast in the form of the 
Helmholtz-Smoluchowski equation U ev = e^Eo/rj which 
describes the electrophoretic velocity for a charged parti- 
cle in the presence of an external electric field [20| , where 
here Eq is given by equation (7). 

We assume the flow is axisymmetric and thus solve 
the system over a two-dimensional cross-section of the 
3-D problem. The 100/im x 100/im simulation domain is 
discretized into approximately 181,000 triangular mesh 
elements. The length and diameter of the simulated 
nanomotor are set to 2 /im and 370 nm, respectively. 
We solve the system of governing equations (1-4) nu- 
merically using the linear system solver PARDISO. Us- 
ing a two-ion model with a fixed bulk electrolyte con- 
centration, the only free parameters in the system are 
the nanomotor native surface potential £ and the sur- 
face flux j. We normalize the flux by characteristic flux 
based on Nernst's diffusion limited current density given 
by j d = 4D+ Coo /\ D [HI. 

Figure 2(a) shows the dimensionless proton concen- 
tration c* = (c+ — Coo)/ Coo (color) and contours of the 
electric potential normalized by the thermal voltage </>* = 
zFcj)/RT (black lines) for the case where ^ = —10 mV 
and j/jd = 0.8. In the absence of reactions, the EDL 
proton concentration and electrical potential are both 
symmetric around the nanorod. Figure 2(a) shows that 
when reaction-driven fluxes are introduced, an asymme- 
try in the proton concentration is established such that 
an excess of protons builds up at the anode and protons 
are depleted at the cathode. The reactions also result in 
an asymmetric electrical potential profile that bulges at 



the cathode. 

Figure 2(b) shows the normalized charge density and 
streamlines of electric field for the same case as Fig. 2(a). 
The charge density in the diffuse layer of the anode is 
positive because the negatively charged surface attracts 
cations near the surface and the surface reactions con- 
stantly inject cations. At the cathode the deficiency of 
protons due to reactions and the shielding protons due 
to the negatively charged surface nearly counteract each 
other resulting in weak negative charge density. The 
charge density generates an electric field, as described 
mathematically by Poisson's equation. The electric field 
couples with the charge density to produce an electri- 
cal body force, p e E, which acts on the fluid to drive an 
electroviscous velocity and propel the nanomotor. 

Figure 2(c) shows the RICA velocity magnitude (color) 
and streamlines (black lines). These simulations are con- 
ducted in the reference frame of the nanomotor. Fluid 
flows from the anode to the cathode due to electrical body 
forces that result from a coupling of positive charge den- 
sity and electric field tangent to the nanomotor surface. 
By Galilean invariance, this is equivalent to the nanomo- 
tor swimming with the anode end forward. 

Figure 3 shows the nanomotor velocity as a function 
of the flux j/jd obtained from simulations, the scaling 
analysis, and experiments. The experiments and sim- 
ulations show good agreement for a native surface po- 
tential of —25 mV, which we obtain from independent 
measurements of the zeta potential for Au and Pt parti- 
cles in aqueous solutions pjj [HI • The nanomotors used 
in the experiments were grown using electrochemical de- 
position The dimensions of these nanomotors were 
similar to those simulated here. We obtained the exper- 
imental data in Fig. 3 by measuring the velocity of ap- 
proximately 200 nanomotors using optical microscopy in 
varying concentrations of hydrogen peroxide. The value 
of j/jd at each peroxide concentration is estimated from 
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FIG. 3: Nanomotor velocity as a function of dimensionless 
flux 2 1 H- Simulations (open symbols), the scaling analy- 
sis (lines), and experiments (closed symbols) show excellent 
agreement. Simulations are shown for four values of the zeta 
potential. 

the published dependence of electrocatalytically gener- 
ated current density on Pt and Au interdigitated micro- 
electrodes [l8| . We have subtracted out the characteristic 
Brownian velocity of the nanomotors (measured here to 
be 4.87 /im s _1 ) from all experimental data points in or- 
der to only consider the axial velocities measured in the 
experiments. 

Additional simulations (not shown) show that the 
nanomotor velocity scales inversely with 77 and D and 



directly with the body force as predicted by the scaling 
analysis. Viscosity slows the nanomotor due to Stokes 
drag and diffusion tends to reduce the proton concen- 
tration gradients generated by the reactions thereby re- 
ducing the generated electric fields. We compute the 
total propulsive force by numerically integrating the z- 
component of electrical body force, p e E z , over the entire 
simulation domain and find that the velocity scales lin- 
early with this body force, in agreement with Stokes law. 
We calculate a force of 0.17 pN at j/j d = 1.0, C = -20 
mV which is in good agreement with our previous exper- 
imental measurement of nanomotor propulsive force of 
0.16 pN [|. 

We have presented detailed simulations, scaling anal- 
ysis, and experiments that describe the locomotion of 
bimetallic nanomotors in hydrogen peroxide solutions 
due to Reaction Induced Charge Auto-Electrophoresis. 
Nanomotor movement is the result of an electroviscous 
slip velocity that is driven by electrical body forces re- 
sulting from charge density and electric fields that are 
internally generated by electrochemical reactions occur- 
ring on the particle surface. We expect that a detailed 
understanding of the physics underlying the nanomotors' 
motion will provide a basis for rational design of next- 
generation nanomachines capable of operation in diverse 
conditions and applications. 

The authors acknowledge Kamil Salloum, Marcus Her- 
rmann, and Joseph Wang for stimulating discussions. 
This work was sponsored by NSF graduate fellowships 
to JLM and PMW and grant CBET-0853379. 



[1] J. R. Howse, R. A. L. Jones, A. J. Ryan, T. Gough, 
R. Vafabakhsh, and R. Golestanian, Phys. Rev. Let. 99 
(2007). 

[2] J. G. Gibbs and Y. P. Zhao, Appl. Phys. Let. 94 (2009). 
[3] S. Fournier-Bidoz, A. C. Arsenault, I. Manners, and 

G. A. Ozin, Chem. Comm. pp. 441-443 (2005). 
[4] J. Burdick, R. Laocharoensuk, P. M. Wheat, J. D. Pos- 

ner, and J. Wang, J. Amer. Chem. Soc. 130, 8164 (2008). 
[5] W. F. Paxton, A. Sen, and T. E. Mallouk, Chem. A Euro. 

J 11, 6462 (2005). 
[6] R. Golestanian, T. B. Liverpool, and A. Ajdari, Phys. 

Rev. Let. 94 (2005). 
[7] N. Mano and A. Heller, J. Amer. Chem. Soc. 127, 11574 

(2005). 

[8] U. M. Cordova- Figueroa and J. F. Brady, Phys. Rev. Let. 
100 (2008). 

[9] Y. Wang, R. M. Hernandez, J. Bartlett, J. M. Bingham, 
T. R. Kline, A. Sen, and T. E. Mallouk, Langmuir 22, 
10451 (2006). 

[10] W. F. Paxton, K. C. Kistler, C. C. Olmeda, A. Sen, 
S. K. S. Angelo, Y. Y. Cao, T. E. Mallouk, P. E. Lam- 
mert, and V. H. Crespi, J. Amer. Chem. Soc. 126, 13424 
(2004). 

[11] R. Golestanian, T. B. Liverpool, and A. Ajdari, New 

Journal of Physics 9, 126 (2007). 
[12] R. Laocharoensuk, J. Burdick, and J. Wang, ACS nano 



2, 1069 (2008). 

T. R. Kline, W. F. Paxton, T. E. Mallouk, and A. Sen, 
Anggewandte Chemie-Inter. 44, 744 (2005). 
P. Calvo-Marzal, K. M. Manesh, D. Kagan, S. Bala- 
subramanian, M. Cardona, G. Flechsig, J. Posner, and 
J. Wang, Chem. Comm. (2009). 

S. Balasubramanian, D. Kagan, K. M. Manesh, P. Calvo- 
Marzal, G. Flechsig, and J. Wang, Small 5, 1569 (2009). 
D. Kagan, P. Calvo-Marzal, S. Balasubramanian, S. Sat- 
tayasamitsathit, K. M. Manesh, G. Flechsig, and 
J. Wang, J. Amer. Chem. Soc. (2009). 
J. F. Hoburg and J. R. Melcher, J. Fluid Mech. 73, 333 
(1976). 

W. F. Paxton, P. T. Baker, T. R. Kline, Y. Wang, T. E. 
Mallouk, and A. Sen, J. Amer. Chem. Soc. 128, 14881 
(2006). 

G. M. Dougherty, K. A. Rose, J. B. H. Tok, S. S. Pannu, 
F. Y. S. Chuang, M. Y. Sha, G. Chakarova, and S. G. 
Penn, Electrophoresis 29, 1131 (2008). 
D. C. Henry, Proc. Royal Soc. Lon. A 133, 106 (1931). 
M. Z. Bazant, K. T. Chu, and B. J. Bayly, SIAM J. App. 
Math. 65, 1463 (2005). 
[22] H. Y. Du, C. H. Wang, H. C. Hsu, S. T. Chang, U. S. 
Chen, S. C. Yen, L. C. Chen, H. C. Shih, and K. H. Chen, 
Diam. Rel. Mat. 17, 535 (2008). 



13 



14 



15 



16 



17" 



18 



19 



